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Abstract 

Electric and color neutral solutions, and the critical conditions for 
the formation of gapless color superconductors, are investigated in K° 
condensed color-flavor locked quark matter for nonzero strange quark 
mass. We show that as the strange quark mass increases, gapless 
modes for up-strange quark pairing occur first, followed by down- 
strange quark pairing. The behavior of the gaps, the dispersion rela- 
tions, and the thermodynamic potential are all found as functions of 
the strange quark mass on the basis of a Nambu Jona-Lasinio type 
model. To a high degree of accuracy, they are presented as relatively 
simple elementary functions. This allows for easy computation for 
any reasonable range of baryon chemical potential and strange quark 
mass. 

PACS numbers: 25.75.Nq, 12.39.Fe, 12.38.-t 
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1 Introduction 



Dense quark matter has been an exciting topic for theoretical study ever 
since it was realized that quark color superconducting phases could form 
the lowest energy states and could be reliably computed at very high den- 
sity within QCD. For three flavors, the original color and flavor S77(3) co ior x 
SU(3) L x SU(3) R symmetries of QCD are broken down to the diagonal sub- 
group SU(3) co \ ot+ l+r at very high baryon density [1]. Quark matter with 
such a particular symmetry pattern is called color-flavor locked (CFL) mat- 
ter and is widely believed to be the densest phase of strongly interacting 
matter [2]. In the situation where the strange quark mass m s is large and/or 
the quark chemical potential ji (which is 1/3 of the baryon chemical poten- 
tial) is small, the so-called gapless color-flavor locked phase (gCFL) has been 
predicted [3, 4, 5]. Similar to the gapless phase of two-flavor color supercon- 
ductor [6] , gCFL is triggered by mismatches between the chemical potentials 
for paired quasi-quarks. There the mismatch in the flavor (say d and s) chem- 
ical potentials is characterized directly by the quantity m^/2/x while that in 
the color (say b and g) chemical potentials is related to m 2 s /2fi via the elec- 
tric/color neutrality in CFL matter [4]. As long as m 2 s /2fi is large enough, 
there exist unpaired (gapless) modes for bd — gs pairing, and the resulting 
gCFL phase becomes energetically favorable. Recently the Meissner masses 
for some gluons in this phase were found to be imaginary so that gCFL is 
actually unstable in the chromomagnetic sense [7, 8, 9]. The presence of a 
chromomagnetic instability is a serious problem. Hence the gapless phase of 
CFL matter needs to be studied further. 

There is another kind of less-symmetric phase of CFL matter that was 
predicted when m s ^ 0. At the leading order, it is convenient to regard the 
effect of m 2 /2/j, as an effective chemical potential associated with strangeness, 
namely 



As one of the pseudo-Goldstone bosons related to excitations of the CFL 
superconducting quark matter (not to be confused with excitations of the 
QCD vacuum with the same quantum numbers), the neutral kaon mode 
K° has a chemical potential /ik° = fJ>s- The mass of this mode, ignoring 
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instanton effects, is [10] 



m° K = — v/3m„(m s + m d ) (2) 

where A is the superconducting gap, m u , rrid and m s are the quark masses, 
and /„- is the pion-decay constant in the CFL matter. The latter is given by 

[11] 

ft = 21 36 ^ 2 ln V - (0.2086^) 2 . (3) 

It measures the strength of the coupling of the pion (in CFL matter) to the 
axial vector current. When ixg reaches the mass of the neutral kaon-mode, 
the K° condenses in CFL matter [10, 12]. This phase is referred to as CFL.fr . 
Assuming rrid <C m s , the condition that no K° condensation takes place is 

m s < 3.034 (m u A 2 ) 1/3 = 44.4 

Taking typical current-quark values of the masses, namely m u = 5 MeV and 
m s = 110 MeV, kaon condensation will occur for reasonable values of the 
gap. In this paper we will generally ignore the up and down quark masses in 
comparison with the strange quark mass. 

Starting with the ideal situation of m s = and raising m 2 s /2ji gradually, 
one expects the CFL matter to be disrupted, at first by the presence of 
K° condensation, and then by the appearance of gapless modes. In the 
CFL matter with K° condensation, the formation of gapless modes has been 
reexamined by using an effective theory in Ref. [13, 14] and by using the 
NJL model in Refs. [15, 16]. It was found that gapless pairing is delayed 
with respect to that in the conventional CFL matter. More recently, it has 
been suggested that the chromomagnetic instability in the gapless phase can 
be resolved by the formation of the so-called kaon supercurrent state [17]. 

Motivated by these results, we investigate the formation of gapless modes 
in the CFLi^ environment. First, we consider the electric/color neutral 
solution of CFLK in the presence of an electron chemical potential and point 
out that the above-mentioned delay is a direct consequence of the deviation 
of the CFLK neutral solution from the CFL one. Second, it is found that 
gapless modes for bu — rs pairing occur first, followed by bd — gs in the K° 
condensed environment. Based on this feature, the resulting gapless phase 
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(termed as gCFLi^ below) is studied, including its electric/color neutrality 
and gap variation. The solutions can be found analytically to very good 
approximation, obviating the need for a complicated numerical minimization 
of the thermodynamic potential with respect to three chemical potentials 
plus three gaps. Finally, the stability of the gCFLfT phase is examined 
qualitatively; it is argued that the gCFL instability might be removed either 
partially or even totally. 

2 Electric/Color Neutrality in Superconduct- 
ing States 

At asymptotically high baryon density the gap in CFL matter can be com- 
puted using weak coupling methods. The result is [18] 




(vr 2 + 4) 



where Nf is the number of flavors. (In this paper we consider Nf = 3 only.) In 
the asymptotic limit, fi ^> Aqcd, it is also true that fi ^> m s or, equivalently, 
^ < A. The only location in the known universe where superconducting 
quark matter might appear is in neutron stars where the baryon chemical 
potential is not asymptotically large. Then jis is not negligible, and further- 
more weak coupling methods in the gauge coupling g are not quantitatively 
reliable. Therefore different groups have used other methods to understand 
what happens in the non-asymptotic regime. These include: (i) variations 
on the Nambu Jona-Lasinio model, which allow a determination of the gaps 
as functions of density at the expense of one or more cut-off parameters, (ii) 
effective Lagrangians employing the relevant collective degrees of freedom 
and with coefficients matched to QCD, but which assume a given value for 
the gap, and (iii) so-called model independent approaches that also assume 
that the gap or gaps are determined by other means and that count states, 
degrees of freedom, and balance electric and color charges. These different 
approaches all tend to agree when they have a region of overlap. 

Stable bulk matter must be electrically neutral and must have zero net 
color. Perfect CFL matter with zero quark masses satisfies these requirements 
without the need for electrons. Alford and Rajagopal [19] addressed this issue 
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in the context of a nonzero strange quark mass by using an essentially model- 
independent approach. They introduced chemical potentials /i e , 113 and /xg for 
the electrons (or negative of the electric charge Q) and for the color generators 
T 3 and 2T 8 / y/3, respectively (The factor of 2/y/3 arises because they used a 
non-standard form of the generator Tg which was more convenient for their 
purposes.) Different flavors and colors of quarks then have different chemical 
potentials; see Table I. The thermodynamic potential Q must satisfy 

dn on on 

«- = «- = «- = °- 6 

0fl e 0/13 Ofls 

They found the solution 

to first order in fi s / fi. (The electron chemical potential is not determined 
by this calculation.) This solution is equivalent to requiring that the quarks 
that should pair have the same "Fermi momentum" to the extent that such a 
quantity has meaning in superconducting matter. To order m\ (and to zero 
order in g 2 ) the difference in thermodynamic potential between the CFL 
phase and ordinary unpaired quark matter is 

Q neutral n ncu tral , ^A* /,,2 a a 2\ /o\ 
CFL - 4 Unpaired + ^ [Vs " 4A ) • ^> 

The CFL phase is preferred over the unpaired phase if the gap is large enough, 
namely if |//£ < A. Another important feature of CFL is that its thermo- 
dynamic potential is independent of 

A*q = -| (l*e + V3 + |a*s) , (9) 

since none of the CFL pairings break the rotation Q = Q — T3 — -h^T%. In 

this sense, CFL matter is not merely an electric insulator [20] but also a Q 
insulator [19]. Enforcing neutrality at higher order requires that no electrons 
be present so that fj, e — 0. 

As mentioned in the introduction, kaon condensation will occur when the 
kaon chemical potential equals its mass. (We must keep in mind that familiar 
names like kaon refer to collective excitations of the color superconducting 
state and not to excitations of the QCD vacuum.) Since we are taking the up 



5 



and down quark masses to be zero (actually only m u = is required in this 
regard) kaon condensation will always occur, and CFL.fr matter should have 
a lower free energy than CFL. In the presence of K° condensation, the ground 
state of Nambu-Goldstone bosons takes non-trivial values. By treating the 
gauge fields as dynamical degrees of freedom, Kryjevski has shown that there 
exist nonzero expectation values for the gluon fields in CFLfT matter [21]. 
He found that electric and color neutrality leads to constant time components 
of two color-diagonal gluon fields, which turn out to be 

°* = t' (10) 

Constant values of the time-components of gauge fields play the mathematical 
role of shifts in the chemical potentials. In this case it corresponds to 

A*3 = -iA*s, V8 = -\ps (11) 

when the condition /i e = is used. (Neutral kaon condensation does not 
require the presence of electrons for electric neutrality.) These results were 
confirmed by Kryjevski and Yamada [14], as were the results of Alford and 
Rajagopal. A comparison of the effective potentials for the quarks in the 
CFK and CFLK states are given in Table I. The thermodynamic potential 
in the CFLi^ state is lower than that of the CFL state (including the strange 
quark mass effects, electric and color neutrality) by an amount \ fl^% [10, 14]. 
Together with (8) this means that 

= nsSSd + % (/4 - 4A 2 ) - \flnl . (12) 

The condition for pairing is now somewhat relaxed to 0.4224/is < A as 
opposed to 0.5^5 < A. 



3 Dispersion Relations 

The modes of excitation in superconducting quark matter are sometimes 
labeled by the quark degrees of freedom and sometimes by the baryon degrees 
of freedom. The correspondence is related to the pairing ansatz 

(q?C>y 5 q$) ~ Aje°*%i + A n e af32 e ij2 + A m e Q/ % 3 • (13) 
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(Here and (a, (3) denote the flavor indices (u, d, s) and the color indices 
(r,g,b), respectively, and the gap parameters A/, Ah and Am are approx- 
imately equal to A in the absence of gap less phenomena.) For example, a 
linear combination of (ru — gd)bu and (gu — rd)bu, where the angular brackets 
denote a pairing, represent an object with baryon number 1, electric charge 
1, and zero net color, in other words a proton. Hence a proton is associated 
with bu quarks. Similarly the neutron, E + , E~, S° and 2~ are associated 
with the bd, gu, rd, gs and rs quarks, respectively. A A , A 8 , or S° would 
be represented by some linear combination of ru, gd and bs quarks. These 
associations are given in Table IV. 

In the CFL phase with massless up, down and strange quarks the three 
gaps are equal: A/ = An = Am = A . The dispersion relations for the 
energy as a function of momentum take the following forms. There are 8 
modes with energy 



The mode with gap 2A is a mixture of ru, gd and bs quarks. 

In the CFL phase with strange quarks of nonzero mass the dispersion 
relations are more complicated because of the necessity of nonzero effective 
chemical potentials that ensure color neutrality. These are found as the zeros 
of a determinant of a 9 x 9 matrix [4] . There are a pair of modes associated 
with each of the pairings bd — gs, bu — rs and gu — rd described by the gaps 
Aj, An and Am, respectively. We shall refer to these as the I, II and III 
channels for brevity. In this phase all three gaps are equal to a common value 
denoted by A, which in general is not equal to A when the strange quark 
mass is nonzero. The pair of dispersion relations for channel I can be written 
as 



involved and fis is one-half of the difference of their chemical potentials; see 
Table II. The same is true in channel II. The pair of dispersion relations for 
channel III are 




(14) 



and 1 mode with energy 




(15) 





(17) 
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The other three modes are linear combinations of ru, gd and bs quarks. There 
are 2 modes with energy 



(18) 



and 1 mode with energy 



E = y/(p - Jif + (2A): 



(19) 



In the CFL.fr phase it is sometimes convenient to use the notation of 
baryons instead of quarks. The translation between the two is given in Table 
IV. The dispersion relations are different from the CFL matter not only 
because of the difference in effective chemical potentials but also because of 
flavor required by the symmetries of QCD [13, 14]. There is a pair of modes 
in channel I with energy given by 



E 



There is a pair of modes in channel II with energy given by 
E = 



2 3 



There is a pair of modes in channel III with energy given by 
E = 



p-(ji+ \ns) 



2 1 



(20) 



(21) 



(22) 



The neutral baryons n, S°, E°, A and A 8 (or bd, gs and ru, gd, bs) mix 
among themselves to produce the following dispersion relations. There is a 
pair of modes with energies 



E= y/(p-fi)* + A*, 



(23) 



and 1 mode with energy 



E = \](p- p) 2 + (2A) 2 . 



(24) 
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4 Transition to a Gapless Phase 



As pointed out in Refs.[4, 6], the gapless phase arises when the energy of 
one of the collective modes becomes zero. As one goes deeper into that 
phase there will be a finite window of momentum for which the energy of 
the mode would be negative. This is the blocking region. One needs to 
find a solution to the combined set of electric/color neutrality conditions and 
gap equations. When gapless modes appear, the gaps for various pairings 
separate from each other and their values should be solved from three gap 
equations, namely dQ/dAj = dQ/dAjj = dQ/dAju = 0. If only one mode 
becomes gapless, it is natural to focus on variations of that gap from the 
original value A while ignoring variations in the other two gaps. If two 
modes become gapless simultaneously then we should consider variations of 
two gaps while ignoring variations of the third. 

The gapless phase will arise when the mismatch between the chemical 
potentials, <5/i, of two quarks exceeds the gap. These mismatches are dis- 
played in Table II for both the CFL and the CFLK phases. See also the 
previous section. The first mode to go gapless in CFL occurs when /i 5 = A 
and is associated with both Aj and Ajj. The first mode to go gapless in 
CFLi^ occurs when fi$ = 4A /3 and is associated with Ajj. We will focus 
on gCFL-fT since we already know that CFLK is favored over CFL when 
(is > 0. 

In Nambu Jona-Lasinio type models the thermodynamic potential is writ- 
ten as 

where Q is a constant adjusted to make the energy of the true vacuum zero. 
The constant G is a coupling constant. The sum runs over 3 flavors, 3 colors, 
and particle plus antiparticle; spin is already included. The energies for the 
particles are generically written as 

Ej (p) = ^(p-^ + A] + % . (26) 

The chemical potentials ftj and 5/ij are linear combinations of //, /is, /is and 
\x e . Anti-particles have chemical potentials opposite to particles. The mo- 
mentum integrals are quartically divergent; usually they are regulated with 



9 



a hard momentum cut-off A. In some sense this approximates asymptotic 
freedom of QCD in the context of this model. The exact expression for the 
integral with such a cut-off is given in the Appendix for completeness. The 
approximate expression which includes terms that go as positive powers of 
A or which remain finite as A — > oo is also given in the appendix. There 
are several terms that are odd in the chemical potentials. These terms can- 
cel when adding particles and anti-particles, so in this sense the inclusion 
of anti-particles simplifies the analysis. Keeping only those terms which are 
finite as A — > oo gives us 



4tt 2 



E 



j-- 



47T 2 127T 2 G 




(27) 

The term proportional to A 4 is absorbed into the vacuum energy. The terms 
quadratic and logarithmic in A are usually accommodated as follows: One 
sets the strange quark mass to zero and solves for the common gap A at 
some reference chemical potential /i for some fixed values of G and A. A 
change in A is accompanied by a corresponding change in G such that the 
gap does not change. Generally it is found that the results are not sensitive 
to the specific choice of A so long as it is large compared to \i. 



4.1 Ideal CFL 

Consider the ideal CFL matter with massless quarks. There are 8 modes 
with gap A and one mode with gap 2 A. These gaps are predicted to be 
very small compared to the chemical potential so the terms in the thermo- 
dynamic potential of order A 4 are normally dropped. The flo and the G get 
renormalized as 

9A 4 

1 A 2 1 

G-l? ~+ G ' (28) 
For convenience we also define 

A' = 2 2 /V 1 A . (29) 
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Then the thermodynamic potential can be written as 

„ ^ V 3A 2 6 2 a2] /A ,s 
127^ + ~G ("A 



(30) 



The solutions to the gap equation, dfl/dA = 0, are either A = or A = Ao 
with the latter determined by 



1 

G 



7T Z 



Ao, 



(31) 



The oft-used requirement that a change in A' should result in a change in G 
such that A be left unchanged means that 

n = 9l (12) 

l + (2//2G /7r 2 )ln(A'/A' ) 1 ) 

So in this case we can trade the pair G,A' for one parameter, namely, Ao- 
Then we obtain the elegant formula 



n = n 



9/i 4 
12^2 



7T^ 



2 In 



A 
Aq~ 



- 1 



(33) 



This obviously has a minimum at A = A . 



4.2 CFL with m s ^ 

Next we consider CFL with fi s ^ 0. The thermodynamic potential now takes 
the form 



n = n 



3 4 , 3A 2 
H 1 



471-2 



7T Z 



2^1„(£)-2^1„(£ 



(34) 



where, as before // = // — The solution to the gap equation is either 

A = or the solution to the equation 

A\ o i /An 



/ , 2 ln(§)=^ln(^)+l(6 / ,- / , 5 K. 



(35) 



When us ^ the dependence on the cut-off no longer drops out. In the limit 
that fis <C n the solution to the above equation is approximately 

A' 



A = An 1 - 



3/x 



2 In 



Ao. 



- 1 



(36) 
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Hence the gap decreases linearly with fis/n- 

Label the solution to the gap equation (35) as Acfl(/^s)- This is a 
definite function of pis- In terms of it the thermodynamic potential may be 
written 

n = sh - % + -Va 2 



2 In | A . ) - 1 

VA C fl(a«5)/ 



(37) 



This obviously has a minimum at A = Acfl(a*s); the dependence on A' is 
implicit through the solution to (35). 

Let us compare with Alford, Kouvaris, and Rajagopal [5] who solved the 
gap equation numerically without approximation. They used Ao = 25 MeV, 
pi = 500 MeV, and A = 800 MeV (A' = 467.18 MeV). In the CFL phase the 
first gapless mode appears when pi s = A. From eq. (36) we find that this 
occurs when A = fi s ~ 23.1 MeV. This is consistent with their Fig. 1. 

4.3 Gapless CFL 

As us increases, eventually one or more of the quark pairings will have a 
mismatch of chemical potentials that becomes as large as the gap. This 
stress breaks apart some of the pairs and reduces that value of the gap. Let 
us write the dispersion relation of quasi-quarks that first become gapless in 
the form 



E(p) = ^(p-fLy + Ai-5fi. (38) 
The excitation energy for quasi-quarks becomes zero at the momenta 

p ± = ^±VV-A 2 . (39) 

This is appropriate when 8 pi > A. Then there exist gapless modes in the 
blocking region p e (p_,p + ). As a consequence, the gapless modes provide 

n,(M*x,A) = J^^[y/(p-fi) 2 + ^-6fA 
5fi 



<5/i 2 - A 2 (5A 2 - 25/i 2 - 12/i 2 ' 

( A -2 A 2\i f S ^+ ~ 



to the total free energy. This must be added to (27). At this point Alford, 
Kouvaris and Rajagopal solve the equations for gCFL numerically using the 
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NJL model. This is the gCFL phase. The three chemical potentials, /i e , /j.3 
and /is evolve as a function of us differently than in the CFL phase. The 
three gaps are no longer equal but evolve such that A/ < A// < A/77 with A7 
dropping the fastest, A77 lagging somewhat behind, and A777 not deviating 
much from A . This is as expected. At the onset of the gapless phase, 
channel III has no chemical potential mismatch and so the associated gap 
A777 should not change much. The channels I and II both have a mismatch 
equal to /is and so their associated gaps should decrease; A/ decreases faster 
than A// because the latter is retarded by the turn-on of the electric chemical 
potential. We shall go no further into the gCFL since the phase with kaon 
condensation has lower free energy. 



4.4 CFLK 

Here we repeat the analysis of section 4.2 but with the chemical potentials and 
dispersion relations appropriate to the CFLK phase. The thermodynamic 
potential is 



1 2 \ 3A 2 



A' 



ln 1 A^ J _ ln 1 A j ~ 



(41) 



where, as before p, = /i — We have dropped terms of fourth order in A 
and us- The functional dependence on the gap is identical to the CFL case. 
Therefore the gap equation and its solutions are also identical. We can thus 
write the thermodynamic potential in the form 



r,2 



3/i 



Art 2 



fi 2 



1 2 



7T" 



2 In 




(42) 



4.5 Gapless CFLK 

Finally we turn our attention to the gCFLA' phase. The difference in free 
energy between the CFL and CFLfT phases was already quoted in (12). We 
now face the problem of how to construct the free energy for the gCFLAT 
phase. The above NJL model allows for different values of the gaps in the 
different channels, and it has the dynamics to determine these in terms of the 
coupling constant G and the momentum cut-off A, but it does not allow for 
kaon condensation. A NJL model was constructed by Forbes [15] which does 
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include kaon condensation. (See also [16].) The resulting calculations involve 
extensive numerical computation, but only the location of the onset of the 
gCFLK phase was studied. The effective field theory that describes kaon 
condensation assumes a common gap and that this gap is determined by other 
considerations. In principle an effective field theory should be constructed 
that makes allowance for differing gaps; this is a formidable task that has yet 
to be worked out. We are interested in describing the gCFLfT phase itself, 
at least for a small range of /is, and not just the location of its onset. 

We shall construct a thermodynamic potential from the following pieces. 
The first piece is the NJL form (25). Practically we will use the large cut- 
off limit (27) and drop the terms of order A 4 . The second piece is the 
contribution of the gapless modes described by (40) which is, however, really 
a part of the NJL model. The third and final piece is the kaon condensation 
free energy. We should allow for the possibility that the values of the chemical 
potentials may be different once the gapless phase is entered. We take this 
from Kryjevski [21]. 



(43) 



Here we use the transcription 



eA°^fi e , gA° 3 ^-fi 3 , gA° 8 



7T 8 - 



This contribution is normalized such that it reflects the free energy difference 
between the CFLK and the CFL phases, namely, 

tt K o(He = 0,/i 3 = -lns,IMi = -\ns) = 

Putting everything together, our model for the CFLi^°/gCFLK matter is 
represented by the function 



3 U 2 111 

O — O ncutral _L 7^2, ,2 _| P a 2 

"CFLifO - ii unpaired + H s + — 

i I 



2 In 



^CFL 



J TV 
III 



~ 1 



|/4 + | (^3 + \PS - He) 2 + | (fJ-8 + \HS ~ §A*c 



+ 'Ee(8tf-/s*)n a (fi i ,5 t i i ,A i )- 
i=i 



12vr 2 



(44) 
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The function AcfiX/^s) is the solution to eq. (35) whose approximate solution 
is given by eq. (36). The average and difference of chemical potentials in 
the i'th channel, fli and 8/ii, are the functions of /i e , /i 3 , fi 8 to be found 
in Tables II and III. The six variables /i e , /i 3 , /i 8 , Aj, A n , and A in must 
be determined on the basis of this function. Note that the solutions give 
the physical thermodynamic potential of eq. (12) when there are no gapless 
modes. This function allows us to study the transition from CFLfT to 
gCFL/T . Note that the gapless mode in channel II is what will first drive 
the chemical potentials and the gap away from their values in the CFLfT 
phase. Before proceeding we will make one other useful approximation; we 
set Jii = p. This has a minimal effect because it is the difference Sfii which 
drives the gapless phase, not the average. In addition we use p, in place of \i 
when evaluating f^. Although f w has not been determined when the quark 
masses are nonzero, this replacement has minimal numerical impact. 

Consider what happens as fis increases, or equivalently what happens as 
the baryon density decreases. Prior to the onset of the gapless modes, we have 
CFL.fr matter. The chemical potentials are /i e = 0,/i3 = — \ns,Hs = — jf^s- 
All three gaps are equal, A/ = A n = A in = A C fl(^s), and decrease with 
increasing ^s- The thermodynamic potential is given by 

= O + % (/4 - 4A^ FL ) - \m . (45) 

The channel II goes gapless when jis = f A CFL (/U,s) ~ 1.206A . The 
difference between the factor of 4/3 (first obtained by Kryjevski and Yamada 
[14]) and 1.206 (very close to the value obtained by Forbes [15]) is due to the 
decrease in the gap with increasing fi S . To study this phase we need to solve 
the equations 

dQ C FLK°/dii 3 = dn CFLK o/dfx 8 = dtt CFLK o/d[i e = (46) 

to obtain the chemical potentials. Allowing for the possibility that any of 
the channels may be gapless, these can be recast in the form 

fl 2/x s - 4/x e + 3^ + 2/i 8 = -— ±- + t;j^ + ~ 2 47 

2 do [in 2 do nin n 2 

^-2^ + 2^1 = ^-^-^ (48) 
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/^ S _2^ + 4 W ] = --^--^ (49) 

Here we have calculated the derivatives db(iijd(i a (where a — e, 3, 8) by using 
Table II. 

When entering the gapless region it is sometimes useful to employ the 
limiting form of Q g when y/5(i' 2 — A 2 <C 5/i, A <^ p,, specifically 

We shall henceforth use this expression. Its validity will be checked a pos- 
teriori. Using this approximation, the gap in channel i must satisfy either 
Aj = or the equation 



h G&)+F ! ^(«-4)-o- 

We are left with solving the above set of six equations to find the behavior 
of (i e , fj,3, fis, A/, An and Am as functions of (is- At this point it is already 
clear that (i e will increase from and a gap will decrease from Acfl (/As) m 
a closely coupled way once we enter a gapless region. 

The electron chemical potential (i e will become nonzero once exceeds 



Acfl- Therefore it is sensible to look for a power series expansion of the 

1/ 



unknowns in terms of x = § (is — Acfl- The results are 




X 3 



/is = -^s + x + w + 



2/i8 = -^ s+x+ h{^ + ^fl) + 



Aj 



^CFL 
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— A CFL ■ 


A m 


= A C FL 




= \»s- 




= AcFL ' 




= A C FL ' 



X 3 



12/i 2 + 



x 3 
x 3 

+ ■ 



12fx 2 
1 x 3 ( 1 



12 U 2 2tt 2 / 2 



+ ■■■ (52) 



The gaps Aj = Am = AcfiX/As) are undeviated from the CFL.fr phase. 
The relative order of magnitude of the corrections displayed above, before 
the ellipses, with x 10 MeV and fx « 500 MeV, is 10~ 4 , indicating that this 
is a very good expansion. This also supports the validity of using the ap- 
proximation (50). To the order displayed in eq. (52) the important quantity 
y/Sfjifi — A 2 H vanishes. It can be determined to lowest nonvanishing order to 
be 

2x 3 



»,-A!, = 5j?+ - (53, 



by inserting the expressions for the chemical potentials from (52) into (48) 
or (49). 

The behavior of /i3, //g and functions of fig are shown in Fig. 1 

while the behavior of Aj = Am and An are shown in Fig. 2. To obtain 
these figures we used eqs. (35) and (52). We used the same numerical values 
for the parameters as in section 4.2, namely, A = 25 MeV, fx = 500 MeV, 
and A = 800 MeV. Essentially all of the functions plotted are linear in fx,s as 
is apparent from the figures. For /is < f Acfl ~ 1.206Ao ~ 30.143 MeV all 
three channels are fully gapped. For larger values of /is channel II becomes 
gapless. 

Channel I becomes gapless when 5/xi = A CFL , which occurs when 

AcFL -^ S "T536V^ + --- 
The chemical potentials and gaps acquire shifts proportional to 

= \ (Us - x 3 /l2^p w f - Aq FL 
V \^ s -x 3 /Wfl)-{^PJ^)A c ^ 1 } 
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Starting from (52) they are shifted as follows: 



a*3 -> - f y + 



2/i 8 -> 2/j 8 + fy + -- • 
vr 2 f 2 

A m -> A m + ... 

tfytj -> - |y H 

-> H 

-> 5jUm + |y H • (55) 

To lowest nonvanishing order 

^tf-Aj = ^y (56) 



Unlike the case of gapless channel II, the corrections here in channel I are not 
suppressed by a factor of I /ft 2 . They begin to deviate rather substantially 
from Acfl- However, we do not need to go very deeply into the gapless 
channel I phase before the color superconducting matter becomes unstable 
to complete breakup into unpaired quark matter. 

Channel I becomes gapless when /j, s > 2A CFL ~ 1.728A xs 43.2 MeV. 
In the region where only channel II is gapless, A n < A/ = A in = A CFL , 
but not by much; it is not visible on the Fig. 2. In the region where both I 
and II are gapless, A/ drops much faster than Ajj. 

With the parameters we have used, channel II becomes gapless when the 
strange quark mass m s = ^/2/i/is > 174 MeV, and channel I becomes gapless 
when m s > 208 MeV. 

The thermodynamic potential, as measured relative to the thermody- 
namic potential of neutral unpaired quark matter, in units of S^Aq/tt 2 , is 
shown in Fig. 3. There is a smooth transition as first channel II and then 
channel I become gapless. Paired quark matter becomes unstable to breakup 
into unpaired quark matter when /is > 50.47 MeV. Coincidentally, this is 
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approximately equal to 2Ao, the value one would get by ignoring the evolu- 
tion of the gaps with us, the chemical potentials enforcing electric and color 
neutrality, and the formation of gapless modes. 

5 Chromomagnetic Instabilities 

Finally, we analyze qualitatively the possibility of eliminating the chromo- 
magnetic instability that occurs in gCFL. In the gCFL case, channels I and 
II become gapless simultaneously at the point (/x e , /is) = (0, Acfl)- The fact 
that the gs — bd and rs — bu pairings become breached simultaneously has 
an important consequence. As stressed in [7, 8], this feature is responsible 
for the instability which occurs for Ai and A 2 gluons. Since the self-energy 
for Ai and A 2 stems from the loop diagram composed of gs and rs quarks, 
the coexistence of gapless (unpaired) gs and rs modes provides a singular 
contribution to the corresponding Meissner masses. In the gCFLfT case, 
however, the gapless formation in channel I is delayed relative to channel 
II. Therefore, singularities should no longer exist in the self-energy for A\^ 
Meissner masses, and the corresponding instability should disappear in the 
gCFLfsT phase. 

On the other hand, the feature that the gCFL location is very close to 
its critical line in the \i e — \is plane has been pointed out to be primarily 
responsible for the instability of A 3 and A 8 gluons and A^ photons [8]. For 
instance, the (almost) infinitesimal channel II gapless strength in gCFL leads 
to the singularities in the A 3 ^ n Meissner masses [8]. In the gCFLK case, 
however, the location is not limited to the vicinity of the critical line. Instead, 
a finite region in the fi e — fi s plane is available for the gCFLif existence. 
Thus such kind of instability might no longer appear, at least for some region 
in the fi e — fi$ plane. 

Since only A 12 and A 3 8)7 exhibit imaginary Meissner masses in gCFL, 
the previous- predicted instability is very likely to be removed in the gCFLK 
phase. Physically, it can be understood from the viewpoint that the chro- 
momagnetic instability is eliminated by building on the proper vacuum. It 
has been suggested that nonzero vacuum expectation values of gluons such 
as (A®) and/or (A$) are helpful for removing the gCFL instability [7]. In 
[22], the instability in two-flavor superconductor was argued to be resolved 
by gluon condensation. In the present work, the nonzero (A®) and (A^) 
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(anisotropic vacuum) derived from the kaon condensation have been at- 
tributed to the change of the electric/color neutral solution and then been 
considered in the study of the gapless formation. In this sense, the conjec- 
ture that gCFL.fr is not chromomagnetically unstable is not very surprising. 
However, intrinsic links between the p-wave kaon supercurrent phase estab- 
lised in [17] and the presently discussed gCFL.fr phase still remain to be 
clarified. Also, we cannot rule out the possibility that other instabilities, 
especially for At, A 5 , A$ and A 7 gluons, arise once again. The quantitative 
calculation of the gluonic Meissner masses in gCFLiT is still necessary, but 
this is beyond the scope of the present work. 

6 Conclusion 

In summary, we investigated electric/color neutrality and gapless formation 
in the CFL matter with K° condensation. By taking the CFLK neutral 
solution into account, we clarify why the gCFLK formation is delayed in 
comparison with the gCFL case. More importantly, it is found that the gap- 
less phenomenon for down-strange pairing (channel I) is delayed while that 
for up-strange pairing (channel II) occurs first. The novel phase structure 
implies that the previously predicted gCFL instability might be removed or 
at least weakened. These conclusions are likely to be important for fully un- 
derstanding the unconventional CFL phases in the presence of strange quarks 
with nonzero masses mj and eletrons fi e . 

Of course, there are still some unanswered aspects in the present work. 
First of all, the K° condensed CFL matter has been treated as a background 
for the gapless formation, which is safe only if /i is not too small; otherwise, 
the condensation may be supressed by instanton effect [23]. The gCFLiT 
phase is lower in free energy than the gCFL phase. Thus, our suggested 
gapless phase should replace the gCFL in neutron star cores [24] , ignoring any 
possible gCFLiT instability. Secondly, more physics involving Goldstone- 
mode condensation should be addressed. These include the possibilities of 
charged kaon and other meson condensations. For instance, if /i e is large 
enough and K~ condensation occurs, the gapless phenomenon for up-down 
pairing in channel III needs to be included. Some of the above-mentioned 
problems are being investigated. 
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Appendix 

An integral that appears frequently is 

jf A dppV (p - PY + A 2 = i/i [6/i 2 + A 2 ] V/i 2 + A 2 



(A -ft) 2 (A - Jif + 4/i 2 + A 2 



+ 3 



(A-/i) 2 + A' 



3/2 



v/(A-/i) 2 + A 2 
^ + A^ 



+^A 2 (4^ 2 - A 2 ) In 



A-^+ J(A-ftY + A 2 



8 v-- y-^ V/U 2 + A 2 - A* 
The limit of this expression as A — > oo is 

/ o A rfp P V(p-f') 2 + A2 - j A " - ^ a3 + j a2a2 + ^ A2A 



2A 



+ A 2 "-/!/ 
+-U (2/i 2 - 13A 2 ) y^ 2 + A 2 



- -/i 2 A 2 + -A 4 
\T T 32 



24' 

The terms which are odd in p, cancel when particles and anti-particles are 
added together. 
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Table I: Effective quark chemical potentials. The third and fourth columns 
evaluate them in the electron-free (fi e = 0) CFL (/i 3 = 0, /ig = —fJ>s) an d 
CFLK (/i 3 = -±Hs,H8 = -\ns) phases. 





effective chemical potential 


CFL 


CFLK 


ru 


H-\n e + \n$ + l(i s 


H- \n s 


li- \li s 


gd 




li- liis 




bs 


\l + \\l e — |/i 8 — \ls 


11 - \lls 


li — f /4s 


rd 


H+\n e + \ns + 5// 8 


li - Ins 




gu 


2 1,1 
A 1 — — 2^ 3 ' 3^ 8 


11- \ii s 




rs 


A* + \Ve + \v>2, + |a*8 - A*S 


11- l/i 5 




bu 


2 2 
A 1 ~~ 3/^2 — 3/^8 


11+ liis 


11+ \ii s 


gs 


H+\He- \n$ + l/l S - fJ-S 


4 

I 1 — 3/^5 




bd 


fl+lUe- \ns 


11+ liis 


/i + i/i 5 
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Table II: Differences of effective quark chemical potentials relevant to pair- 
ing, defined as 5/i — \ (/i q i — /i^) for the pair ql-q2. The third and fourth 
columns evaluate them in the electron-free (/i e = 0) CFL (/i 3 = 0, a 8 = —u s ) 
and CFLK (u 3 = -|/x s ,// 8 = —\fJ>s) phases. 





5u 


CFL 


CFLK 


ru — gd 


\ (-a e + /i 3 ) 







ru — bs 


\ ("/ie + \^3 + ^8 + a S ) 





\»s 


gd — bs 


\ (-lV3 + H8 + Vs) 





\l*s 


bd — gs 


\ (5^3 - A*8 + A*s) 


As 


\l*s 


bu — rs 


| (-/i e - 5^3 - A*8 + /4s) 


Ms 


\»s 


gu — rd 


I ("/ie - /i 3 ) 








Table III: Averages of effective quark chemical potentials relevant to pairing. 
The third and fourth columns evaluate them in the electron-free (u e = 0) 
CFL (u 3 = 0,u 8 = -u s ) and CFLK (u 3 = -|/x 5 ,/x 8 = -\fis) phases. 





average 


CFL 


CFLK 


ru — gd 




H-\ns 


A* - J2^S 


ru — bs 


V - \He + \[J>3 - \V-8 - \l*S 


H-\ns 


A* - J2^S 


gd — bs 


U + |/i e — |/i 3 — — |/i5 


a - lu s 


a - l(j, s 


bd — gs 


U + |/i e — \a 3 — |/i 8 — \u s 


A 4 — I Us 


H — §A4s 


bu — rs 


V - \ll>e + \^3 - |a*8 - \HS 


H-\ns 


A* - i2^S 


gu — rd 


V ~ \He + \^8 


a - \a s 
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Table IV: Correspondence between quark and baryon descriptions with 
condensate form (qfC^q?) ~ Aje Q/ %-i + A H e a ^e ij2 + A ni e al33 e ij3 . 





baryon label 


ru 


mixture of E°, A , A 8 


gd 


mixture of E°, A , A 8 


bs 


mixture of E°, A , A 8 


gu 


E+ 


rs 




bu 


P 


gs 


wo 


bd 


n 
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Figure 1: Chemical potentials as functions of /is- The numerical parameters 
chosen are // = 500 MeV, A = 25 MeV, and A = 800 MeV. The curves are 
terminated when the paired phase becomes unstable to decay to the unpaired 
phase. 
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Figure 2: The gaps as functions of us- The numerical parameters chosen are 
Ii = 500 MeV, A = 25 MeV, and A = 800 MeV. The curves are terminated 
when the paired phase becomes unstable to decay to the unpaired phase. 
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Figure 3: The difference in thermodynamic potentials between the paired 
and unpaired phases as a function of us- The numerical parameters chosen 
are fi = 500 MeV, A = 25 MeV, and A = 800 MeV. The transition occurs 
when /is ~ 50.47 MeV. 
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